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Abstract 

A  general  and  qualitatively  exact  theory  is  developed  for  quantum  sticking 
coefficients  <x(k)  in  the  small  wave  number  limit  k  -►  0.  The  theory  covers 
Morse-type  to  inverse-square  potentials,  the  latter  representing  long-range 
potentials.  The  theory  gives  unambiguous  answers  to  crucial  questions  in  the 
problem  and  helps  lead  to  an  overall  understanding  of  low- temperature 

adsorption.  f  ~Tc^.  :  > :  L(ir 
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I .  Introduction 

The  problem  of  determining  sticking  coefficients  of  an  atom  on  material 

surfaces  is  one  of  the  fundamental  tasks  in  surface  physics.  A  general 

situation  is  depicted  in  Fig.  la.  If  the  kinetic  energy  of  incident  atoms  or 

molecules  is  large  compared  with  the  energy  scale  of  the  problem,  e.g.,  a 

potential-well  depth  of  matter-particle  interaction,  then  it  is  widely  believed 

that  the  distorted-wave-Born-approximation  (DWBA)  is  quantitatively  correct  for 
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describing  various  surface  phenomena.  ’  On  the  other  hand,  when  the  incident 
kinetic  energy  is  much  smaller  than  the  typical  energy  scale  of  the  given 
problem,  as  realized  in  low  temperatures,  there  has  been  a  long-standing 
controversy  concerning  the  validity  of  the  DWBA.^  ^  This  is  partly  due  to  the 
lack  of  accurate  experiments  at  small  k  and  at  low  temperatures.  The  only 

g 

exception  is  the  experiment  done  by  the  Edward  group. 

The  present  theoretical  situation  is  the  following.  People  believe  that  a 
one-dimensional  (1-d)  simplification,  as  depicted  in  Fig.  lb,  is  good  enough  to 
study  qualitative  properties  of  the  sticking  coefficients  a(k)  at  low 
temperatures.  It  is  noted  that  a  T-shape  model  as  shown  in  Fig.  lc  has  a 
deficiency  that  there  is  no  momentum  sink  along  the  direction  of  the  particle 
motion.  In  classical  mechanics,  a  particle  would  eventually  strike  the  surface 
no  matter  how  slow  it  is,  releasing  some  energy  to  a  matter  system,  and 
therefore,  it  would  definitely  be  adsorbed  on  the  surface  for  vanishingly  small 
incident  energy.  Quantum  mechanically,  however,  the  situation  is  not  quite  so 
obvious.  First,  it  appears  natural  that  in  quantum  mechanics  a  low-energy  wave 
function  of  the  atom  must  have  a  small  wave  amplitude  at  a  high  potential 
barrier,  i.e.,  at  a  surface.  This  is  true  for  static  short-range  potentials  due 
to  the  following  reason.  Consider  a  static  square-well  potential  as  shown  in 
Fig.  2.  The  scattering  wave  function  with  incident  wave  number  k  can  be  written 
as 


for  0  <  x  <  R 
for  R  <  x  , 


(1.1) 
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where  S(k)  is  an  S-matrix  element  and 


K2K2  h2k2 
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where  m  is  the  particle  mass.  Imposing  the  smoothness  condition  on  |k+>  at 
x  *  R,  one  can  easily  find  that 


(1.2) 


A  «  k  ,  S(k)  -*•  1  for  k  0 


(1.3) 
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The  linear  k-dependence  of  the  scattering  wave  amplitude  in  the  potential-well 
region  immediately  leads  to  a  linear  k-dependence  of  the  sticking  coefficient 
a(k)  at  small  k.  The  first  quantum-mechanical  calculation  of  a(k)  by  Lennard- 
Jones  and  Devonshire  (UD)  based  on  the  DWBA  is  essentially  the  same  as  the 
above  argument  for  the  static  square-well  potential.  Based  on  the  LJD 
calculation  for  a  static  Morse  potential  which  decays  exponentially  at  large 
distance,  people  next  considered  various  possible  mechanisms  for  realizing 
finite  a(0).  These  include  static  long-range  potentials,  x  and  correlated 
or  many-body  motions  of  the  atom-matter  system  such  as  polaron^  and  self- 
trapping,^  but  for  short-range  potentials  like  square-well  and  Morse-type. 
However,  all  the  existing  theories  share  a  common  deficiency  in  that  they  do  not 
answer  the  following  two  questions:  How  important  is  the  long-range  character 
of  the  matter-particle  interaction?  How  important  is  the  many-body  motion  of 
the  matter-particle  system?  It  is  clear  that  a  satisfactory  theory  must  treat 
these  effects  on  an  equal-footing.  Are  there  any  potentials  for  which  a(0)  are 
finite?  If  there  are  such  potentials,  what  is  a  borderline  potential  which 
distinguishes  between  "short-range"  potentials  for  which  a(0)  ■  0  and  "long- 
range"  potentials  for  which  a(0)  *  0? 


.VsVsVaVv-  ■ 
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In  this  paper,  we  shall  develop  a  new  theory  of  quantum  sticking 

coefficients  at  0  K.  We  do  not  invoke  any  uncontrolled  approximations,  and  thus 

give  a  definitive  answer  to  the  problem.  Our  answer  is  as  follows:  (i) 

Irrespective  of  the  matter-particle  interaction  potentials,  the  many-body  motion 

of  matter-particle  system  is  not  essential  for  the  small-k  behavior  of  a(k). 

The  present  theory  gives  a  firm  basis  to  the  DWBA.  (ii)  For  potentials 

_2 

decaying  faster  than  x  at  large  distance,  a(k)  «  k  at  small  k.  This  is 

essentially  due  to  the  fact  that  the  scattering  wave  function  has  a  linear  k- 

dependence  at  small  k.  The  range  of  the  linear  k  region  decreases  with 

increasing  long-range  character  of  the  matter-particle  interactions.  The  theory 

•  - 1  4 

encourages  experiments  at  small  k  <  0.1  A  .  For  a  He  atom  colliding  with  a 
liquid  He  surface,  the  present  theory  predicts  a  critical  wave  number  -  0.01 
A  *  below  which  o(k)  depends  on  k  approximately  linearly,  which  provides  a  deep 
understanding  of  the  experiment.  (iii)  The  inverse-square  potential  is  a 
borderline  distinguishing  between  short-range  potentials  for  whish  a(0)  •  0  and 
long-range  potentials  for  which  a(0)  *  0.  This  conclusion  is  reached  by  a  DWBA 
calculation  of  a(k)  for  the  inverse  square  potential.  Our  assertion  (i)  in  this 
case  means  that  the  DWBA  result  cannot  be  changed  qualitatively  by  higher-order 
corrections.  This  is  essentially  due  to  the  fact  that,  unlike  the  case  of 
short-range  potentials,  the  scattering  wave  function  here  has  a  /E-singularity 
at  small  k. 

Before  getting  into  detailed  descriptions,  let  us  try  an  intuitive 

explanation  of  the  present  theory  and  its  main  results.  To  do  so,  consider  once 

more  the  static  square-well  potential.  Fig.  2.  In  the  above,  we  have  given  a 

standard  explanation  of  a(0)  ■  0.  Our  method  in  this  simple  case  would  proceed 

as  follows.  We  first  note  that  the  S-matrix  element  S(k)  ■  exp(2i6(k)),  where 

q 

4(k)  is  a  scattering  phase  shift.  But  Levinson's  theorem  tells  us  that  6(k)  is 
related  to  the  number  of  bound  state  N  of  the  square-well  potential  as 
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4(+0)  -  «N  .  (1.4) 

Sow  for  k  <<  1,  S(k)  -  1  by  (1.4)  and  R  -  0(1)  for  short-range  potentials,  and 
therefore 

|k+>  a  -  2ikx  for  R  <  x  <<  l/k  .  (1.5) 

On  the  other  hand,  |k+>  in  the  region  0  <  x  <  R  must  be  a  linear  combination  of 

two  independent  solutions  of  the  Schrodinger  equation,  e**^.  Noting  that  K 

does  not  depend  on  k  in  the  limit  k  0,  and  requiring  a  smooth  connection  of 

|k+>  at  x  *  R,  one  can  easily  see  that  both  two  constant  factors  before  e**^ 

are  proportional  to  k,  which  leads  to  the  same  result  as  before  --  |k+>  «  k  at 

the  potential  well,  and  hence  a(k)  «  k  for  small  k.  Now  the  point  is,  this 

explanation  does  not  use  the  boundary  condition  at  x  ■  0,  nor  does  it  require 

details  of  the  wave  function  in  the  potential  well.  Therefore,  under  the 

existence  of  a  dynamical  version  of  Levinson's  theorem,  it  can  be  extended  to 

dynamical  cases,  which  is  precisely  what  we  ^hall  do  in  the  next  section.  The 

existence  of  a  dynamical  version  of  Levinson's  theorem  for  matter-particle 

-2 

interaction  potentials  decaying  faster  than  x  has  recently  been  shown  by  the 
present  authors.^ 

We  have  organized  the  present  paper  as  follows.  In  the  next  section, 
short-range  potentials  are  examined  based  on  formal  scattering  theory  and  the 
dynamical  Levinson's  theorem.  In  Section  III,  a  numerical  calculation  of  a(k) 
is  carried  out  based  on  a  self -trapping  model,  which  demonstrates  the  dynamical 
Levinson's  theorem  and  determines  a  standard  number  p  to  be  introduced  in 
Section  II.  In  Section  IV,  we  examine  the  inverse  square  potential  which  is 
essentially  different  from  the  short-range  cases  and  is  outside  the  scope  of  the 
arguments  to  be  presented  in  Section  II.  Our  DWBA  calculation  presents,  for  the 
first  time,  an  example  of  finite  a(0).  Our  summary  is  given  in  Section  V. 
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(After  the  completion  of  the  present  work,  we  have  noticed  the  paper  by  Brivio  ! 

\ 

and  Brimley  in  which  the  adsorption  problem  is  treated  by  the  Green's  function 

method.  ^*)  j 

( 

I 

I 

i 

II.  Short-Range  Potentials 

We  first  examine  short-range  potentials,  that  is,  those  decaying  faster  j 

-2 

than  x  .  Let  us  consider  the  adatom  +  1-d  matter  Hamiltonian  in  a  collinear 

I 

configuration  (see  Fig.  lb), 

Hfcot  *  H(x,p)  +  V(x,x)  +  K(p)  ,  (2.1) 

where  H(x,p)  is  a  1-d  matter  Hamiltonian  with  x  and  p,  respectively, 

representing  position  and  momentum  vectors  for  atoms  in  the  1-d  matter,  V(x,x) 

is  an  interaction  potential  between  the  adatom  and  1-d  matter  which  we  restrict 

-2 

to  those  decaying  faster  than  x  at  large  adatom-matter  separation  x  >>  1,  and 
K(p)  is  the  kinetic  energy  of  adatom.  It  is  noted  that  the  1-d  matter  must  be 
of  finite  size,  because  otherwise  only  one  phase  is  possible  in  one  dimension, 
that  is,  there  is  no  distinction  between  gases,  liquids  and  solids.**  This 
means  that  the  matter  does  not  have  a  well-defined  boundary  at  finite 
temperatures,  and  the  question  of  calculating  an  adsorption  probability  becomes 
meaningless. 

A  scattering  eigenstate  of  (2.1)  with  a  wave  number  k  can  be  written  as 

♦(?,k)  -  F(r,k)  -  S(k)F(r , -k)  ,  (2.2) 

••  •  •  “  ilex  «• 

where  r  s  (x,x),  F(r,k)  is  a  Jost  solution  having  an  asymptotic  form  e  $q(x) 

it  x  +  *,  where  ^q(x)  is  the  ground  state  of  H(x,p)  (T  *  0  K),  and  S(k)  ■ 

exp[2i6(k))  is  an  S-matrix  element.  In  a  recent  paper, ***  we  have  extended 

Levinson's  theorem  in  static  potential  scattering  to  our  dynamical  case.  We 


have  shown  that  for  our  short-range  potential  V(x,x),  the  scattering  phase  shift 
6(k)  is  connected  to  the  number  of  bound  states,  N,  of  the  total 
Hamiltonian  by 
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6(+0)  -  Nit 


(2.3) 


On  the  other  hand,  an  exact  expression  of  an  adsorption  rate  R(k),  the 

number  of  adsorped  atoms  per  unit  time,  is  given  by  time-dependent  scattering 
.v  12 

theory  as 

R(k)  *  Hi  \  6(Ek  ‘  Ef)|<f lVf lk+>|2  »  (2*4) 

f 


where  E,^  = 

(T  -  0  K), 

*  + 
interaction,  e.g.,  creation  of  phonons,  and  |k  >  represents  a  scattering 

eigenstate  of  the  total  Hamiltonian,  Hfcot.  The  sticking  coefficient  a(k)  is 

defined  as  a  ratio  between  R(k)  and  the  flux,  the  number  of  incident  atoms  per 

unit  time,  which  is  given  by  Hk/m.  Therefore  we  have 


h2k2 

+  Eq  with  Eq  representing  the  ground  state  energy  of  the  matter 
T  means  a  summation  over  final  states  |f>,  is  a  final  channel 


a(k) 


(2.5) 


An  important  observation  in  (2.5)  is  that  as  far  as  the  small-k  behavior  of  a(k) 
is  concerned,  it  is  given  essentially  by  the  quantity  |k  >  / k.  This  is  seen  as 
follows.  The  static  potential 


V(xq,x)  *  V(x,x)  -  Vj  ,  (2.6) 

where  Xq  denotes  a  mean  configuration  of  matter  atoms  in  the  ground  state  $q(x), 
is  by  assumption  a  short-range  potential,  and  supports  only  a  finite  number  of 
bound  states  with  finite  binding  energies.  Moreover,  the  final  states  |f>  in 
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(2.5)  are  eigenstates  of  the  unperturbed  Hamiltonian  H  -V^.  Therefore,  in  the 
limit  k  -+  0,  all  the  quantities  other  than  |k  >  /k  in  (2.5)  do  not  depend  on  k. 

We  now  explore  the  small-k  behavior  of  the  exact  scattering  eigenstate 
|k+>.  From  the  dynamical  Levinson's  theorem  (2.3),  for  k  <<  1  the  S-matrix 
element  S(k)  =  exp[2i6(k)]  -  1  in  (2.2).  Let  us  define  a  reduced  potential 


V(xQ,x)/(^)  =  U(x)  . 


(2.7) 


and  consider  large  distances  R  and  R'  such  that 


-U(R)  *  k«  -  U(R')  «  { -U(  x) } 


(2.8) 


This  means  that  the  matter  system  is  still  in  the  ground  state  in  the  region 

**  *  ilex  — 

R'  <  x,  and  that  the  Jost  solution  in  (2.2)  is  F(r,k)  s  e  $q(x)  t^ie  region 
R  <  x.  Thus  we  have 


1,+^  ,  * ikx  ikx. 

|k  >  s  (e  -  e  )  <*q(x) 


(2.9) 


s  -2ikx$g(x)  for  R  <  x  <  -  , 


where  in  the  second  expression,  we  have  assumed  the  inequality 


kR  <  1. 


(2.10) 


Next  consider  the  region  R'  <  x  <  R,  where  the  particle  motion  is  described  by 
Schrodinger  equation 


{-  +  U(x)}  *  •  k2*  . 

dx 


(2.11) 


Instead  of  explicitly  solving  (2.11)  for  respective  potentials  U(x),  however,  we 
simply  point  out  that  the  particle  motion  in  this  region  must  be  somewhere 
between  the  two  extreme  cases,  U(x)  *  U(R)  and  U(x)  =  U(R'),  as  shown  by  two 
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horizontal  dashed  lines  in  Fig.  3.  But  one  can  easily  see  that  in  both  cases 
the  smoothness  requirement  of  |k+>  at  x  =  R  [cf.  (2.9)]  leads  to  |k+>  «  k  at  R' 
<  x  <  R.  This  fact  in  turn  leads  to  |k+>  “  k  at  0  <  x  <  R',  because  in  this 
region,  |k+>  can  depend  on  k  only  through  a  constant  prefactor  which,  however, 
again  by  smoothness  of  |k+>  at  x  =  R',  must  be  proportional  to  k.  In  this  way 
we  reach 

|k+>  «  k  0  <  x  <  R.  (2.12) 

From  (2.12)  we  have 

a(k)  «  |k+>2/k  «  k  .  (2.13) 

We  remember,  however,  that  (2.13)  is  derived  under  the  condition  (2.10)  and 
S(k)-  1.  Little  is  known  about  the  small-k  dependence  of  S(k).  Our  analysis 

below  and  numerical  calculations  in  the  next  section  for  a  Morse  potential  show 

that  the  latter  condition  imposes  a  similar  restriction  on  k  as  (2.10).  To 
analyze  (2.10)  for  typical  long-distance  behaviors  of  the  reduced  potential 
U(x),  let  us  introduce  a  standard  number  3  <  1  in  analogy  to  the  "standard 
value"  in  the  Lindemann  formula  for  melting,*2  and  replace  (2.10)  by 

kR  <  3  .  (2.14) 

With  this  3  we  can  predict  for  each  given  potential  a  critical  wave  number  k£ 
such  that 

a(k)  «  k  at  k  <  k  .  (2.15) 

-  c 

2 

From  (2.14)  and  U(R)  ■  k  ,  one  can  easily  obtain  the  following  results: 

(i)  Exponential  decay  U(x)  -  -ye  *X  , 
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-(k  /ic)ln(k  /y)  =  P 
c  c 


(2.16) 


(ii)  Algebraic  decay  U(x)  -  -yx  n  (n  >  2)  , 


1/n  .  1-2/n  0 

Y  kc  =  8  . 


(2.17) 


The  results  (2.16)  and  (2.17)  make  two  predictions.  First,  the  critical 


wave  number  k£  decreases  with  increasing  long-range  character  of  the  potential. 


Our  numerical  calculations  of  a(k)  in  the  next  section  for  a  Morse  potential 


(2.16)  with  y  -  48  A  *  and  x  *  1.4  A  ^  (numbers  for  the  W-He  system^’^)  show 


that  k  -  0.1  A  .  Substituting  these  numbers  into  (2.16)  gives 


P  -  0.6  . 


(2.18) 


4  4 

The  He  atom  colliding  on  a  liquid  He  surface  is  described  by  a  van  der  Waals 


potential  (2.17)  with  n  *  3  and  y  =  20  A.  Using  (5  -  0.6  in  (2.17),  we  have  k 


*  1  8 
0.01  A  for  this  case,  which  is  in  good  agreement  with  experiment.  Second, 


from  (2.17)  kc  approaches  zero  quite  rapidly  as  the  exponent  n  approaches  2  from 


above,  as  depicted  in  Fig.  4  for  y  =  20  and  (5  =  0.6.  This  strongly  suggests 


finite  a(0)  for  an  inverse  square  potential,  and  that  n  *  2  is  a  borderline 


distinguishing  between  "short-range”  potentials  n  >  2  and  "long-range" 


potentials  n  <  2.  We  note  that  the  n  =  2  specialty  is  beyond  physical 


intuition,  but  is  clear  mathematically  as  signaled  in  several  ways:  the 


dynamical  Levinson's  theorem  (2.3)  breaks  down  for  n  <  2,  the  inequality  (2.10) 


does  not  hold  for  n  <  2,  and  the  zero-energy  scattering  wave  function  has  a 


singularity  at  n  ■  2.  This  last  fact  is  seen  as  follows.  For  k  •  0  and  U(x) 


- Y*  >  (2.11)  becomes 


*"  +  Yx’n<i  »  0 


(2.19) 


A  general  solution  of  (2.19)  is 


v  v  >  /■. 


□ 


(2.20 


* 


A  z  .  (£2  xl-»/2 

1  2-n 

2-n 


where  is  a  solution  of  Bessel's  differential  equation 


Z"  +  -  Z'  +  (1 

V  XV 


0  . 


(2.21) 


The  singularity  of  *  at  n  -  2  is  clear  in  (2.20).  In  Section  IV,  we  shall 
present  a  thorough  study  of  the  inverse  square  potential  and  show  that  the 
quantum  sticking  coefficient  a(k)  is  finite  in  the  limit  k  0,  thus  confirming 
the  above  prediction  that  n  ■  2  is  a  borderline. 


Ill .  Determination  of  B 

In  this  section,  we  consider  a  self -trapping  model  for  the  adsorption  of 
He  atom  on  W.  A  collinear  configuration  of  He  and  W  is  shown  in  Fig.  5.  By 
calculating  the  sticking  coefficient  a(k>  as  a  function  of  the  wave  number  k,  w t 
shall  demonstrate  the  dynamical  Levinson's  theorem  (2.3)  and  determine  the 
standard  number  p.  The  model  assumes  a  harmonic  lattice,  H  *  a*a  a 

Ph  [  qqq 

q 

surface  W  atom  plus  He  atom  interacting  through  a  Morse  potential,  and  the 
surface  W  atom  seeing  another  Morse  potential  created  by  the  end  W  atom  of  the 
1-d  lattice  in  its  equilibrium  (taken  to  be  x  *  0) 

Hs+a  "  Ps/2M  +  0(xs)  +  P»/2“  +  U(xa‘xs)  *  (3.1) 

The  dynamical  interaction  between  the  harmonic  lattice  and  the  surface  W  atom 
takes  the  form  (one-phonon  process) 


H 


int 


h 

M/.. 


)i  (a!  +  a 


d0(xs) 


I 


e(t)  , 


(3.2) 


where  is  the  number  of  lattice  atoms,  and  the  unit  step  function  0(t)  means 
that  |k+>  here  is  a  scattering  eigenstate  of  Hs+a>  not  that  of  the  total 
Hamiltonian,  In  analogy  to  the  second-order  optical  process,  this  is  a 
luminescence  component  approximation  thereby  neglecting  Raman  components,  which 


is  expected  to  be  good  for  short-range  potentials  of  the  Morse  type.  To 
calculate  a(k)  by  (2.5),  we  note  that  here  is  at  t  >  0,  and  the  final 

state  is  |f>  =  a^|0>|b>,  where  |0>  is  the  phonon  vacuum  and  |b>  represents  a 
bound  state  of  H  .  . 


To  solve  the  Schrodinger  equation 


H  di 
s+a 


0.3) 


which  is  the  main  task  in  this  section,  we  now  specify  the  W-He  and  W-W  Morse 
potentials  as 


0(x+x°) 

s 

U( x+x°) 

l  a  J 


(3. A) 


where  x^  -  3  A,  x^  -  3.6  A,  D  -  0.99eV,  D  -  10  meV  and  <  z  1.4  A  The 

s  a 

phonon  dispersion  is  given  in  terms  of  the  W-W  interaction  potential  0  as 

Mw  ^  *  20  (x*x^)[l  -  cos(qx^)]  .  (3.5) 

q  xx  s  ^  s 


It  is  worth  noting  that  the  problem  (3.3)  is  essentially  a  three-body  one  and  an 
analytic  solution  is  not  available  even  for  the  simplest  interaction  potentials, 
e.g.,  the  square-well.  The  best  we  can  do  is  to  construct  perturbative  and 
variational  wave  functions  and  solve  (3.3)  numerically.  In  doing  so,  we  note 
the  following.  First,  the  one-body  problem  for  Morse  potentials  is  exactly 
solvable.  Appendix  A  summarizes  some  results  about  Morse  potentials.  In  the 
case  of  the  W-He  Morse  potential  U(x),  the  quantity  K  as  defined  in  (A. 3)  is 
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about  7.4,  which  by  (A. 6)  means  that  there  are  four  bound  states  with  the  lowest 

-21 

binding  energy  -  -1.52  *  10  J  s  -9  meV.  Similarly  for  the  W-W  Morse 
potential  0(x),  there  are  about  200  bound  states.  Measuring  the  energy  from 
that  of  the  lowest  bound  state,  the  first  four  bound  states  have  energies  0  (n 
*  0),  1.35  (n  =  1),  2.69  (n  =  2)  and  3.94  (n  =  3)  in  the  unit  of  10  ^  J.  The 
one-particle  energy  spectra  for  the  W-W  and  W-He  Morse  potentials  are 
schematically  summarized  in  Fig.  6.  It  is  now  clear  from  energy  considerations 
that  as  long  as  k  is  small,  the  first  three  bound  states  are  good  enough  to 
describe  the  motion  of  the  surface  W-atom.  As  for  the  motion  of  He,  there  are 
two  situations.  When  the  He  atom  is  far  away  from  the  surface,  the  surface  W- 
atom  is  in  the  ground  state  and  the  wave  function  is  written  as 

-ikx  ikx 

*  *  *asym(VXa)  *  C*0(xS)le  3  '  S(k)e  &]  *  (3'6) 

where  ♦Q(xs)  is  the  ground  state  of  the  surface  w-atom  as  given  by  (A. 2),  and 
C  *  0  for  bound  states  |b>  and  C  ■  1  for  the  scattering  state  |k+>.  Near  the 
surface,  on  the  other  hand,  the  wave  function  can  be  written  as 

M-l 

*  '  *surf<V*,)  '  l  8p(X a)V,»>  •  ”'7) 

p*0 

where  the  g^'s  describe  a  fully  correlated  motion  of  He  near  the  surface. 

According  to  the  energy  consideration  above,  we  take  M  *  3  hereafter.  One  may 

suggest  the  use  of  the  variable  x  -x  in  place  of  x  in  (3.7)  and  expand  the  He 

motion  in  terms  of  the  eigenstates  in  (A. 2)  for  the  W-He  Morse  potential 

U.  This  scheme  is  not  efficient,  however,  because  of  complexities  arising  from 

2 

the  kinetic  energy  term  ps/2M  in  (3.1).  We  find  that  the  Gauss ian-weighted 
Hermite  polynomials  are  the  best  suited  orthonormal  complete  set  over  which  to 
expand  localized  g  -functions: 


14 


-T2x2/4 

h  -  8  e  '  H  (Tx) 
q  q  q 


q  ~  0,1,2...  , 


(3.8) 


where  8  ■  [T/ql/Ti)]1  and  the  H  's  are  Hermite  polynomials.  In  Appendix  B, 

q  q 


we  have  described  Hermite  polynomials  and  related  integrals  which  appear  in  the 


calculation  of  the  Hamiltonian  matrix  elements.  The  scale  factor  T  in  (3.8) 


must  be  best  fitted  according  to  the  degree  of  locality  of  the  g^'s.  To  check 


the  efficiency  of  the  h^  expansion,  we  have  expanded  the  first  four  bound  states 


of  the  W-He  Morse  potential  in  terms  of  the  h  's.  We  have  found  that  14  h  's 

q  q 


reproduce  these  four  wave  functions  very  well,  with  <  1Z  relative  error  in  wave 


amplitudes. 


Based  on  these  preliminary  arguments,  we  now  consider  the  following 


variational  wave  functions  for  the  eigenstates  of  Hs+a: 


2  N-l 

<|>(x  ,x  )  *  *  ♦(*.)  +  y  y  c  ♦  (x  )h  (x  )  , 

s  d  asym  a  L  L  ij  i  s  j  a 

i=0  j-0 


(3.9) 


where  *(*a)  is  an  error  function  preventing  the  He-atom  from  penetrating  into 


the  matter. 


rx-x^-x°-/2/v  22 

v  (  as  , .  -v  t 

n  J  dt  e 

J  -00 


(3.10) 


where  v  is  a  variational  parameter.  First  consider  the  bound  states  |b>. 


Substituting  (3.6)  with  C  ■  0  and  (3.9)  into  (3.3)  and  operating 


lbs*'.  WW  • 


we  obtain  a  set  of  coupled  algebraic  equations  for  the  coefficients  C^'s: 


l  V 


C  b  -  ec  b 
ijCij  eClm  * 


(3.11) 


\  *.  V  \  *  A  .>  A.‘»  A  -  ^  .  •  .  •  Aj,  *  -  .*-X  .'j  Va  • 

*_*  ‘JL.  J  A. 


where 


H.  .  .  s  c.60 .6  .  -  6.  .P  .  +  e  0.  .  .  , 

i.m,ij  x  ii  mj  £.1  mj  g  £m,ij 


(3.12) 


and  where  e.  =  2mE./J^,  e  =  2mE/)^,  e  =  2Dm/h^,  and  P  .  and  0  , 

ii  g  mj  fcm.ij 

respectively,  are  given  by  (B.7)  and  (B.9).  Starting  with  (3.6)  and  C  *  1,  one 

can  obtain  a  similar  equation  as  (3.11)  for  the  scattering  state  |k+>.  It  is 

noted  that  unlike  the  case  of  bound  states,  the  eigenenergy  is  already  known  as 
2  2 

=  HI  k  /2m.  However,  there  is  an  additional  unknown,  S(k),  and  some 
additional  matrix  elements  associated  with  this  variable.  We  have  obtained 


H.  .  -  e.6OJ6  .  ' 

£m,ij  k  ti  mj  i 

i 

-H* 

£m 

H*.  ! 

ij  i 

i 

H 

ss 

(3.13) 


where  e,  =  k  ,  and  H  ,  H  and  H  are  certain  integrals  involving  the  error 
lc  Jtni  s  ss 

function  »,  (3.10),  and  are  given  in  Appendix  C. 

We  can  now  calculate  the  sticking  coefficient  a(k)  as  defined  by  (2.5).  As 

noted  before,  our  self -trapping  model  assumes  the  one-phonon  final  states  |f>  * 

a*|0>|b>  and  takes  as  |k+>  the  scattering  eigenstate  of  ^s+a»  and  here  is 

H.  (3.2),  at  t  >  0.  Let 
int 

‘  2ii  b£2  ’  1  “  1’2”,"t0 

be  the  binding  energy  of  the  t-th  bound  state.  With  the  phonon  spectrum  (3.5), 
the  energy  conservation  in  (2.5)  becomes 


5(E*-V  ■  zt  TvR  Ht"«>  +  s<q+q°,)  • 


(3.14) 


where 


(3.15) 


v/M/m  /mB  ()c  +  b^),  qQ  1  -q  cos  1  ( 1  -B^ )  . 


As  for  the  matrix  element  in  (2.5),  from  (3.2),  (3.6)  and  (3.9)  we  have 


KflvflkV  -  fopiV2 

0  q 


(3.16) 


with 


Mt  -  Y  Cy  Jdx  hj*(e'ikx-  S(k)eikX} 


Jdx  ...Oo- 


(x)/2tcC 


ij 


+  l  Cy  l  CJj  Jdx  Vp  °'(x)/2kC 

ij  P 


N  N 

£  (Y(  j  ,0)  -  S(k)Y*(j,0)}-^{KJ(K;i,0;l)  -  J(K;i,0;2)} 

ic  1C 


ij 


N  N 

+  y  Cii  y  CDi  -L“2<KJ(K; i,p;l)  -  J(K; i ,p;2)>  , 

L  J  L  PJ 


ij  P 


(3.17) 


where  (A. 11)  and  (C.6)  are  used.  Substituting  (3.14)  and  (3.16)  into  (2.5) 
gives 


«<k>  -  ^(g)*^)3'2  y - |h  i: 

#  A  ,,  „2  kz  +  b.  1 


(3.18) 


1=1  /2-BT  "  ’  “l 

t 


We 


have  numerically  solved  (3.11)  and  (3.13)  and  obtained  b^,  C^,  S(k)  and 


Cj*.  We  have  then  evaluated  a(k)  based  on  (3.15),  (3.17)  and  (3.18).  The 


stability  and  convergence  of  the  numerical  calculation  has  been  checked  by 
changing  variational  parameters  T  and  v,  as  well  as  the  number  of  Gaussian- 


weighted  Hermite  polynomials,  h^.  To  test  the  validity  of  our  numerical 


calculations,  we  have  first  checked  if  the  numerically  obtained  binding  energies 


vt'jp .  v  m  tt  *.■  v  r  ■.»  ■.■  v  ^ 1  w  w.vjwmw  ^  ^  ^  ^ 

Ri 


£> 

I 
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v 

« . 
l  . 


are  reasonable.  We  have  found  four  bound  states  with  energies  -20.087,  -9.053, 
-2.939  and  -1.837  in  the  unit  of  A  .  Simple  estimates  of  these  energies  are 
the  four  binding  energies  of  the  W-He  Morse  potential,  assuming  the  surface  W- 
atom  to  be  in  the  ground  state  $q.  From  (A. 6),  in  the  same  unit  of  A  ,  these 
estimates  are  -20.133,  -9.530,  -2.846  and  -0.082,  which  are  close  to  the 
numerical  results,  thus  demonstrating  the  validity  of  our  numerical 
calculations.  Figure  7  contains  our  final  results  for  T  *  2,  v  *  1  and  N  *  21. 
Part  a  shows  that  the  S-matrix  element  takes  the  limit  S(k)  ->  1  as  k  +  0, 
demonstrating  the  dynamical  Levinson's  theorem,  (2.3).  Part  b  shows  a(k)  as  a 
function  of  the  wave  number  k.  The  approximate  linear  k-dependence  of  a(k)  is 
seen  below  k  -  0.1  A  Substituting  this  into  (2.16)  with  <  *  1.4  A  *  and  y  - 
48  A  ^  gives  p  -  0.6. 

IV.  Inverse-Square  Potential 

In  the  preceding  two  sections,  we  have  investigated  short-range  cases, 

-2 

i.e.,  the  particle-matter  interaction  potentials  decaying  faster  than  x  at 
large  distance,  x  >>  1 .  Due  to  the  existence  of  the  dynamical  Levinson's 
theorem  (2.3),  we  could  develop  a  general  and  accurate  theory  for  the  small-k 
behavior  of  the  sticking  coefficients  a(k).  As  we  have  found,  a(k)  «  k,  thus 
verifying  that  the  DWBA  is  qualitatively  correct  in  this  case.  One  important 
result  is  that  a(k)  «  k  only  at  k  <  kc>  and  the  critical  wave  number  kc 
decreases  with  increasing  long-range  character  of  the  potentials.  In 
particular,  as  Fig.  4  shows,  the  theory  predicts  that  kc  ■  0  and  hence  a(0)  is 
finite  for  an  inverse-square  potential.  In  other  words,  it  predicts  that  the 
inverse-square  potential  is  a  borderline  distinguishing  between  short-range 
potentials  for  which  a(0)  ■  0  and  long-range  potentials  for  which  a(0)  is 
finite.  In  this  section,  we  shall  carry  out  a  one-phonon  mediated  DWBA 
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calculation  of  the  sticking  coefficient  for  an  inverse  square  potential,  and 
show  that  this  is  indeed  the  case. 

Before  getting  into  details,  it  is  useful  to  look  at  the  problem  of 
sticking  coefficients  from  the  viewpoint  of  formal  perturbation  theory.  The 
exact  scattering  eigenstate  |k+>  in  (2.S)  can  be  written  as 

oo 

|k+>  «  y  (G0Vf)n|k>  ,  (4.1) 

n*0 

where  is  a  dynamical  interaction  between  the  adatom  and  matter,  |k>  is  an 

exact  scattering  eigenstate  for  a  static  adatom-matter  potential,  and  Gq  =  (Ek 

-H  ♦  i0+)  *  is  an  adatom  propagator  with  the  Hamiltonian  H  describing  the  adatom 

moving  in  the  static  adatom-matter  potential.  Substituting  (4.1)  into  (2.5) 

provides  a  power  series  expansion  of  a(k)  with  respect  to  V^,  whose  lowest 

order,  0(V^  ),  is  the  DWBA.  In  the  case  of  short-range  potentials,  the  DWBA 

gives  a  vanishing  a(0),  and  it  has  been  a  challenging  question  if  the  inclusion 

of  higher-order  terms  can  bring  about  a  qualitatively  different  answer,  namely  a 

finite  o(0).  In  the  case  of  the  inverse-square  potential,  on  the  other  hand,  as 

we  shall  see  below,  the  DWBA  already  gives  a  finite  a(0).  Now  an  important 

2 

point  is  that  the  DWBA  is  the  only  term  to  0(V^  ).  This  means  that  an  exact 
a(0)  is  also  finite  in  general.  In  other  words,  as  far  as  the  qualitative 
small-k-behavior  of  a(k)  is  concerned,  there  is  no  significant  many-body  effects 
in  the  case  of  the  inverse-square  potential. 

Another  important  remark  is  on  the  dimensionality  of  the  problem.  As  we 
have  noted  in  Section  I,  a  1-d  model  of  infinite  size  is  physically  meaningless 
for  the  study  of  adsorption.  Therefore,  when  one  uses  the  1-d  model,  the  system 
size  must  be  finite,  so  that  the  phonon  spectrum  has  a  gap.  In  the  case  of 
short-range  potentials,  the  number  of  bound  states  if  finite,  and  therefore  a 


Vv.V  sy. 


cutoff  is  effectively  introduced  in  the  pi.onon  spectrum.  This  is  why  the 
commonly  used  infinite  1-d  model  does  not  show  any  difficulty.  For  the  inverse- 
square  potential,  on  the  other  hand,  as  we  shall  see  below,  there  are  infinitely 
many  shallow  bound  states  which  strongly  couple  to  very  "soft"  phonons,  leading 
to  an  infrared  catastrophe  to  be  explained.  We  thus  encounter  a  finite-size 
problem  in  a  1-d  model  for  the  inverse-square  potential.  To  avoid  this  problem, 
we  shall  extend  the  commonly  used  1-d  model  to  higher  dimensions  --  two  and 
three  --  thereby  expecting  that  the  phonon  density  of  states  at  higher 
dimensions  would  suppress  the  contribution  of  very  soft  phonons  and  remove  the 
infrared  catastrophe. 

Let  us  consider  a  particle  of  mass  m  striking  the  material  surface  in  the 
normal  direction.  The  Hamiltonian  is  written  in  general  as 

N0 

H  *  p2/2m  +  Hph  +  ^  V(r-l)  .  (4.2) 

i=*l 

Writing  iT  *  iTg  +  6fT  and  following  the  standard  procedures,  we  divide  the 
third  term  into  a  static  potential 

v(z)  -  ]v(M*.0)  (4.3) 

i 


and  a  dynamical  interaction 


-  l  1N0(^).?  ei<q^>  r 
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In  (4.3),  we  have  assumed  a  translational  symmetry  parallel  to  the  surface.  In 

(4.4),  e  is  a  phonon  polarization,  V  is  a  Fourier  ii$- component,  and  $  is  a 
qX  <$ 

reciprocal  lattice  vector.  Now  the  one-phonon  mediated  DWBA  for  a  sticking 

coefficient  a(k^  means  that  in  (2.5),  the  final  state  is  Jf>  ■ 
t  ik|fr||  + 

a  |0>|b>e  11  ",  and  |k  >  and  |b>  are,  respectively,  scattering  and  bound 

q\ 

eigenstates  of  the  static  potential  v(z).  Here  the  subscript  ||  denotes  a 
component  parallel  to  the  surface.  From  (2.5)  and  (4.4)  we  have 


-X  l  ^ 


qg  +  b2)  -  #u)  }  |M| 2 

"  qA 


(4.5) 


with 


2N0M», 


r1No  l 

qX  2  *  (0,0, G .) 


i(qz  +  Gz)z  + 
<b|e  |k  > 


x  (q  +  S)-C  V 

qA  q-Hj 


(4.6) 


In  deriving  (4.5)  and  (4.6),  we  have  put  Gj|  *  0  because  at  low  energy  a  nonzero 

Gy  does  not  satisfy  energy  conservation.  Moreover,  in  a  standard  situation, 

such  as  a  He-atom  colliding  on  W-matter,  we  see  from  energy  conservation  that 

2 

q  <<  1,  which  enables  us  to  neglect  the  qii  -term  compared  with  in  (4.5),  and 

•*  qA 

approximate  (4.6)  as 


dsF  "o  <bi  l 
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r*  iG  Z 

)  e  2  iG  t*  V 
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(4.7) 


We  note  that  the  G  *  0  term  in  (4.6)  is  negligibly  small.  On  the  other 

z 

hand,  one  can  easily  show  that 
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v(z)  =  No  l 


iG  Z 
e  2  V,. 


(4.8) 


With  (4.8),  (4.7)  can  be  written  as 


"  *  <>i  <J» db-  ^  'k+> 

i^V  o 


(4.9) 


Here  it  is  noteworthy  that  the  result  (4.9)  justifies  replacing  the  dynamical 
interaction  (4.4)  by  an  effective  one. 


2N-Mu> 
0  ■* 


(a-  +  a  ->  > 

qX  -qX 


dv(  z) 
dz 


(4.10) 


which  is  a  natural  extension  of  the  commonly  used  1-d  model  of  low- temperature 
adrorption  to  higher  dimensions. 

To  explicitly  evaluate  (4.5)  and  (4.9)  for  a(k),  we  must  specify  the  static 

potential  v(z)  and  the  phonon  dispersion  w  .  For  the  static  potential  we  take 

qX 

the  form 


r  -2 
-5  YZ 


2  ^  2o 

Z  <  z„ 


(4.11) 


where  z  *  Zq  denotes  the  surface.  It  is  noted  that  our  phonon  argument  above  is 
not  consistent  with  the  long-range  potential,  (4.11).  Conduction  electrons  will 
be  more  relevant  than  phonons  for  the  long-range  interaction.^  However,  as 
will  become  clear  below,  the  specific  form  of  the  interaction  term  is  not 
essential  for  finiteness  of  a(0),  and  (4.10)  is  good  enough  to  study  qualitative 
aspects  of  quantum  adsorption  for  the  long-range  potential.  For  a  similar 
reason,  we  consider  W  and  He  as  an  example  of  matter  and  atom.  Clearly  the  W-He 
interaction  is  not  of  the  type  (4.11),  and  the  choice  of  W  and  He  here  is  simply 
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for  the  purpose  of  specifying  the  atom  and  the  phonon  dispersion.  We  replace 
u  by  ui  of  (3.5),  which  upon  substitution  of  (3.4)  becomes 


«  -  2<(j-j{l  -  cos(qx^)}]1^2 


(4.12) 


This  is  based  on  the  observation  that  (4.12)  gives  a  sound  velocity  -  4  x 
3 

10  m/s,  which  is  close  to  the  geometrical  mean  of  the  longitudinal  and 

_  3 

transverse  sound  velocities  in  three  dimensions,  -  3.4  '  10  m/s.  With  this 
simplification  and 


[  <*!  )2  ■  1  • 

LX  qA 


(4.13) 


the  dimensionality,  D,  appears  only  in  the  phonon  density  of  states 


D 

?^T  J0  d"  °d(,!) 


(4.14a) 


PD(q)  »  2  (D  «  1),  x°q  (D  -  2)  and  (x°q)2/it  (D  =■  3)  .  (4.14b) 

Our  final  task  towards  the  calculation  of  a(k)  is  to  solve  the  Schrodinger 


equation 


(4.15) 


which  is  carried  out  in  Appendix  D.  Figure  8  contains  the  long-range  potential 
(4.11)  in  the  unit  of  A  2  with  2q  ■  0.5  A  and  p  *  /y  *  l  m  5.0,  and  its  bound 
states  calculated  from  (D.9).  From  (3.14),  (4.5),  (4.9)  and  (4. 1 1 ) - ( 4 . 14),  we 


have 


ihf* 

4W5m, 


I 


=  -X-1 — 2  l<^N'3|k+>|2  . 
k  +  b* 


(4.16) 


where  is  the  i-th  b  as  given  by  (D.9),  and  | fc>  is  the  t-th  bound  state 
(D.14),  with  the  first  bound  state  being  the  lowest  bound  state.  Substituting 
(D.4),  (D. 6) ,  (D.14)  and  (D. 18)  into  (4.16)  gives 


a(k)  = 


4ir)i  y  "irjj 


►VBm  zi 


lJv^z0^  \  PD^q0^ 


,  .2  lI-v+i^btzo)  *  Iv-i(bezo)l  lMJ 


✓2-B?  k  +bt 
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where  v  *  Vi*Y  and 


(4.17) 


dz  z  (J  (kz)J  (kzn)  -  J  (kz)J  (kzn)} 
v  -v  0  -v  v  0 


*  {Iv(biz)  *  I-v(V)}  • 


(4.18) 


For  given  values  of  parameters  Zq,  p  and  k,  one  can  calculate  the  prefactor  in 
(4.17)  and  bound-state  spectrum  b^  by  (D.9),  (D.ll)  and  (D.12),  which  in  turn 

_2 

determines  and  PD(qQ)  by  (3.15)  and  (4.14b),  and  ^^(b^Zg)  ’  *v-l^biZ0^ 
by  (D.ll).  Finally,  the  matrix  element  (4.18)  can  be  evaluated  numerically, 
with  some  remarks  given  in  Appendix  E. 

We  can  now  discuss  some  properties  of  a(k).  First  of  all,  the  finiteness 
of  a(0)  can  be  seen  as  follows.  Consider  the  lowest  bound-state  contribution  in 
(4.16).  It  is  clear  in  this  expression  that  in  the  limit  k  -*■  0,  the  k- 
dependence  of  a(k)  is  determined  by  that  of  the  quantity  |k  >  / k,  which, 
however,  approaches  a  constant  due  to  the  /E-dependence  of  |k+>  in  the  limit  k  •+ 
0,  (D.7).  Therefore,  the  lowest-bound-state  contribution  to  a(0)  and  hence  the 
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value  of  a(0)  itself  is  constant.  An  important  point  here  is  that  this 

conclusion  of  finite  a(0)  does  not  use  details  of  the  interaction  H.  „ .  Thus, 

int 

we  have  seen  that  the  finiteness  of  a(0)  is  a  general  consequence  of  the 
inverse-square  potential.  Secondly,  we  have  observed  a  logarithmic  divergence 
of  a(k)  in  the  1-d  case, 

aCk)^  ,  «  -  Ink  for  k  -+  0  .  (4.19) 

D—  1 

To  see  the  origin  of  this  catastrophe,  consider  the  matrix  element  M  in  the 
limit  k,b  -+  0.  In  this  limit,  substituting  (D.2)  and  (D.ll)  into  (4.18)  gives 


".-C 


dz  z"2Uv  - 


(«.*)  1 


r(l-v)  r -V  1 _  r  -  V  1 

f.  /  i  \  i 


r(l+v)  t  l-2v  l  l+2v  ’ 


(4.20) 


where  5  b^z^/2.  From  (4.20)  and  (D.12)  we  have 


9  9  A 

| M  |  «  -  v  /|l+2v|  =  constant. 


(4.21) 


On  the  other  hand,  in  the  same  limit  k,b  ■*  0,  (D.2),  (3.15),  (D.18)  and  (D.19) 


give 


-*■  constant,  -*■  0 


I’-v+dV  -  Iv-l(Bt)l  *  6t 


-1 


(4.22) 


Substituting  (4.21)  and  (4.22)  into  (4.17)  and  using  the  density  of  bound  states 
(D.13),  we  find  the  contribution  of  the  bound  states  0  <  b  <  bp  <<  1  to  a(k)  as 


°  db  p  (k2  +  b2)  r 
0  U  k^  +  b 


oc 


(4.23) 
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which,  for  D  *  1,  diverges  logarithmically  in  the  limit  k  -*  0.  That  is,  in  the 
limit  k  -*  0,  each  of  the  contributions  from  very  shallow  bound  states  remains 
finite,  and  a(k)  diverges  logarithmically  due  to  the  linearly  diverging  density 
of  bound  states.  However,  as  we  have  pointed  out  in  Section  I,  the  obtained 
infrared  catastrophe  does  not  imply  a  new  surface  many-body  problem,  but  reveals 
a  pathological  aspect  of  the  infinite  1-d  model. 

Finally,  we  have  numerically  evaluated  (4.17)  and  (4.18)  for  the  D  *  1,  2 
and  3  cases.  Figure  9  contains  sticking  coefficients  a(k)  as  functions  of  the 
wave  number  k.  The  parameters  used  are  Zq  =*  0.5  A  and  p  *  5.0.  Note  that  a 
logarithmic  divergence  of  a(k)  in  the  D  *  1  case  is  spurious.  Although  the 
infinite  2-d  model  has  in  general  the  same  pathological  aspect  as  the  1-d  model 
in  that  there  can  be  no  lattice  formations  in  one  and  two  dimensions  just  as 
there  can  be  no  spontaneous  magnetizations  in  the  isotropic  Heisenberg  model  in 
these  dimensions, ^  the  spurious  contribution  of  very  soft  phonons  is  suppressed 
by  the  decreasing  phonon  density  of  states  at  q  -♦  0  [cf.  (4.14b)],  as  a  tiny 
difference  between  a(0)  and  a(2)  in  the  D  *  2  curve  demonstrates.  The  results 
for  the  cases  D  =  2  and  3  demonstrate  finite  a(0)  for  the  inverse-square 
potential . 

V.  Summary 

In  this  paper  we  have  developed  a  new  theory  of  quantum  sticking 

coefficients  at  0  K.  Based  on  the  dynamical  Levinson's  theorem,  we  have  first 

considered  the  cases  where  the  adatom-matter  interaction  potentials  decay  faster 

_2 

than  x  .  We  have  demonstrated  that  the  many-body  motion  of  adatom-matter 
system  is  not  essential  in  determining  the  low-energy  behavior  of  the  sticking 
coefficients.  The  essence  is  instead  in  the  long-distant  part  of  the 
interaction.  The  theory  predicts  that  a(k)  «  k  at  k  <  k£,  and  that  the  critical 
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wave  number  k^  decrease  with  increasing  long-range  character  of  the  adatom- 

matter  interaction.  To  give  a  quantitative  measure  of  kc ,  we  have  introduced 

the  standard  number  p  in  analogy  to  the  "standard  value"  in  the  Lindeman  melting 

formula.  To  determine  p,  we  have  considered  a  self -trapping  model  for  the  W-He 

system.  By  carrying  out  precise  numerical  evaluations  of  a(k),  we  have 

determined  that  p  -  0.6.  As  a  test  of  the  present  theory,  we  have  then  applied 

A  A 

our  results  to  the  case  of  He-atom  colliding  with  a  liquid  He  surface,  for 
which  accurate  experimental  data  are  available.  The  theory  predicts  that  a(k)  « 
k  at  k  <  0.01  A  *,  in  good  agreement  with  experiments.  Next,  we  have  considered 
the  inverse-square  potential.  Our  study  based  on  the  dynamical  Levinson's 
theorem  has  already  verified  that  k^  -*■  0  as  the  potential  approaches  the 
inverse-square  one,  which  strongly  indicates  a  finite  a(0)  for  the  inverse- 
square  potential.  We  have  carried  out  the  one-phonon-mediated  DWBA  for  the 
quantum  sticking  coefficient  at  Q  K  and  have  shown  that  the  inverse- square 
potential  is  a  borderline  distinguishing  between  "short-range"  potentials  for 
which  a(0)  *  0  and  long-range  potentials  for  which  a(0)  *  0.  We  note  that  the 
finite  a(0)  result  of  DWBA  would  hardly  be  changed  qualitatively  by  including 
higher-order  corrections.  It  is  also  noted  that  this  result  depends  on  the  fact 
that  the  scattering  wave  function  has  a  /Ii - s ingular ity  in  the  small-k  limit,  in 
sharp  contrast  to  the  linear  k-dependence  in  the  case  of  short-range  potentials, 
and  therefore  it  is  a  general  consequence  of  the  inverse-square  potential.  In 
brief,  the  essence  of  quantum  sticking  coefficients  at  0  K  is  in  the  k- 
dependence  of  the  scattering  wave  function  |k+>  in  the  low-energy  limit  k  0. 

For  short-range  potentials,  |k+  has  a  linear  k-dependence  which  leads  to  a 
linearly  vanishing  a(k).  On  the  other  hand,  for  long-range  potentials,  |k+>  has 
a  i/Tc-singularity  leading  to  a  finite  a(0). 


Acknowledgments 


This  research  was  supported  by  the  Air  Force  Office  of  Scientific  Research 
(AFSC),  United  States  Air  Force,  under  Contract  F49620-86-C-0009,  the  Office  of 
Naval  Research  and  the  National  Science  Foundation  under  Grant  CHE-8620274.  Thi 
United  States  Government  is  authorized  to  reproduce  and  distribute  reprints 
notwithstanding  any  copyright  notation  hereon. 


28 


Appendix  A.  Some  Results  about  Morse  Potentials 
The  ore-body  problem  for  a  Morse  potential. 


.2 

{‘  2m  2  +  D(e*'<X  '  2e"'°C>>  ♦(*)  =  E*(x)  , 

dx^ 


(A.  1 ) 
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is  exactly  solvable.  The  bound  states  are 
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(A. 2) 


where 


K  =  2/2mD/htc  ,  z  =  Ke 


-  KTX 


(A. 3) 


and  the  F^s  are  Laguerre  polynomials. 
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satisfying  the  differential  equation 
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The  corresponding  binding  energies  are 
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The  normalization  constant  N  is  determined  by  the  orthonormal  condition 
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The  integral  in  (A. 7)  is  a  special  case  (£  *  0)  of  the  Nieto  and  Simmons 


integral 
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where  T  is  a  gamma  function.  In  particular, 

J(K;i,i;0)  -  (K'J_1)  r(K-2i-l)  , 

from  which  the  normalization  constant  is  determined  as 
Nt  *  /Ki ! (K-2T-1 )/(k-i-l) !  . 

Similarly, 

r  KlN  N 

j  dx  KX  =  — - — 1  J(K;i,j;-t)  . 

Appendix  B.  Hermite  Polynomials  and  Related  Integrals 
The  Hermite  polynomials 
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Using  (B.2)-(B.4),  we  obtain 
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Next,  using  the  expansion  (B.l)  we  have 


if"  dx  h  c 
I  m 

*  -00 


-Ax 


Q  (A)  =1  dx  h  e  h 
mn  I  m  n 


A  jj  [m/2]  [n/2] 

£  £  (-l)r+r'(2r-l)l!  (2r'-l)!l 

r*Q  r'*0 


*  r 


.  -x  /2-Ax/T  m-2r+n-2r' 

dx  e  x 


(-A/T)”*"  eA2/2T2 


Cm/2]  [n/2]  [(m+n)/2]-r-r' 


l  l  l 


x  (2r-l ) ! !  (2r'-l)!!  (2k-l)l !  (A/T) 


(-1) 


r+r’ 


r«0  r 1  *0  k“0 

-2r-2r'-2k 


/O  \  /  n  wO+ti-2r-2r' .. 

l2rM2r,M  2k  ' 


(B.5) 


(B.6) 


(B.7) 


X 


(B.8) 


31 


Using  (A. 11)  and  (B.8)  we  evaluate  the  Integral 
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Appendix  C.  Some  Integrals  Involving  ♦ 
First  consider  the  integral 
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With  (3.8),  (3.10)  and  (B.3),  A  becomes 
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where  we  have  defined 
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Using  (3.8)  and  (B.l)  and  carrying  out  a  standard  Gaussian  integral,  we  obtain 
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As  for  the  term  B,  from  (3.4)  and  (A. 11)  it  becomes 
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where  we  have  defined 
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Integrating  by  parts  and  using  (3.8),  (3.10),  (B.3),  and  (B.4),  we  have 
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From  (C.3)  and  (C.7)  we  have  the  recursion  relation 
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The  first  term  Y(0,c)  cannot  be  evaluated  analytically,  but  the  dimension  of 
integrations  can  be  reduced  from  2  to  1.  In  the  expression 


W(c,T,k,v)  3  Y(0,c)  -  J  dx  hQ  ♦  e’(c+lk)x 


Rnv  -|~(x+v^2/u)2  -  (c  +  ik)(x+/2/u)  (x  2  2 

77*  I  dx  e  dt  e  ,  (C.9) 

J  -QB  »  -OD 


we  make  the  variable  change  t  +  x+t,  and  the  order  of  integrations  as 


Carrying  out  a  Gaussian  integral  over  x,  we  reach  the  results 
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Next  consider  integrals  which  involve  the  error  function  ♦  twice 
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s  a  0  as 


2  2  2 
e  KN  2e  KN 

-*  ^  °  J(K;0 ,0 ; -2)  G(k,2<)  +  ^  -°-  J(Kj0,0;-l)G(k,x) 


where 

—  / .  \  -  f  j  *2  (2ik-c)x 

G(k,c)  =  dx  ♦  e 

With  the  use  of  (3.10),  an  integration  by  parts  and  (C.9),  we  have 

G(k,c)  -  c-'2ik  W(c-2/Zv,2v,-2k,v)  • 

We  must  also  carry  out  the  following  integrations: 

H2  =  J  dx  ♦(♦"-2ik*')  *  -  (ik+v//2?) 

u2  -  f  j  *  ’s  2ikx  v  -k2/2v2  +  2/Zik/v 

Hss  =  J  dx  ♦  (*  +2ik*  >  e  "  *  * 

and  Hss  in  (3.13)  are  then  given  by 
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H  =  Hl  (k)  +  H2 
ss  ss  ss 


(C . 1 7 ) 


Appendix  D.  Solution  of  (4.15) 

The  solution  of  (4.15)  can  be  written  as 


♦  *  /z  [AJ^(ez)  +  BJ_^(ez)]  , 


(D.l) 


2  2 

where  y  =  1/A  -v  =  1/4  +  p  ,  and  A  and  B  are  constants  of  integration.  is 

18 

the  Bessel  function  of  the  first  kind, 


J  (z)  *  (z/2)V  V  (-1)  (z/2) _ 

V2'  -  U/21  £  n,r(v  +  n  +  1}  • 


(D.2) 


Since  (4.15)  is  real,  we  can  take  as  real.  First  consider  a  scattering  state, 
where  e  *  k  is  positive.  Noting  that  J*(z)  *  J  ^(z),  we  can  satisfy  the  reality 
condition  and  the  boundary  condition  ^(Zq)  *  0  by  taking 


B*  =  A  *  iJ_v(kzQ)  . 


(D.  3) 


With  (D.l)  and  (D.3)  we  now  write  the  scattering  state  as 


|k+>  »  iN^/z( J^(kz)J_  j(kzn)  -  J_u(kz)Ju(kzn) ) 


■v  0 


v'""'  v'O' 


(D.4) 


To  find  the  normalization  factor  N^,  note  the  asymptotic  property  at  |z|  -*•*»/ 


V2)  -/ScosU  *  *>  • 


(D.5) 


By  setting  the  prefactor  of  the  e  component  in  |k  >  at  z  ■+  «  equal  to  1, 


is  given  by 


Nk  -  /Z5I  |J.v(kzQ)  e'wy/2  -  Jv(kzQ)  e*y/2|"1 


-  ✓25E  e_1TM/2/|Jv(kz0)| 


(D.6) 
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where  in  the  second  expression  we  have  assumed  p  >>  1 .  The  case  p  <  1 
corresponds  to  an  unrealistically  narrow  and  deep  potential.  A  particular 
interest  is  in  the  limit  k  ■+  0,  where  from  (D.2),  (D.4)  and  (D.6)  we  have  for  z 
«  k'1 

|k+>k+Q  =  i/kz/p  [z/zQ)V  -  (z/zQ)  V]  *  A  .  (D. 7 ) 

This  A-singularity  of  the  scattering  state  |k+>  is  in  sharp  contrast  to  the 
linear  k-dependence  of  Jk+>  for  the  short-range  potentials,  (2.12).  As  we  see 
in  the  text,  this  is  the  key  point  of  finite  a(0)  for  the  inverse-square 
potential. 

Next  we  consider  bound  states.  In  this  case,  e  *  ib  (b  *  positive).  The 
reality  condition  is  met  by  taking 

B  -  A*e_1TM  .  (D. 8) 

On  the  other  hand,  the  finiteness  requirement  for  the  wave  function  at  z  +  « 
leads  to  a  pure  imaginary  A.  With  (D.8)  and  the  fact  that  A  is  pure  imaginary, 
the  boundary  condition  iJi(zq)  *  0  provides  the  quantization  of  bound  states, 

W  ‘  I-v(b*0)  *  0  ■  <D9) 
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where  Iy(z)  is  the  modified  Bessel  function  of  the  first  kind. 

Alternatively,  one  can  look  at  the  Jost  function  f(-k),  the  prefactor  of  the 
•  component  in  the  scattering  state  |k+>  at  : 

f(-k)  -  Jv(kzQ)  e1^72  -  J_v(kzQ)  e_iry/2  ,  (D.10) 

whose  zeroes  at  k  ■  ib  in  the  complex  k-plane  describe  the  bound  states.  The 
condition  f(-ib)  *  0  is  nothing  but  (D.9).  Unlike  the  short-range  potentials 
studied  in  Section  II,  the  inverse-square  potential  has  an  infinite  number  of 


WWW 


v.v  . 


»V  ij.i » ft  »_y«  7*  r* '.'»  rnrerr  rv  vvr  w  yv 
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bound  states.  To  see  this,  consider  the  limit  bz^  -♦  0.  In  this  limit,  only  the 

1 8 

n  *  0  term  in  the  power  series  expansion  of  I^(z), 


Iv(z) 


(,/2)''  v _ uw2n 

(zl2)  l  n! r( v  +  n  +  1)  * 

n=0 


is  kept  and  the  quantization  condition  (D.9)  becomes 


JL  lnlLl±yl  .in 

2  2v  1  T(l-v)  p 
—  e  e 

zo 


(D. 11) 


(D. 12) 


where  n  is  any  integer  such  that  bZg  <<  1.  That  is,  there  are  inf initely-many 
shallow  bound-states,  whose  density  of  states  diverges  at  b  0  as 

p(b)  *  p/nb  .  (D. 13) 

Finally,  to  normalize  the  bound  states,  we  use  (D.8)  and  ReA  *  0  to  write 

| b>  *  iNfe/z  (Iv(bz)  -  I^(bz)]  =  iNb/E  [-  2sin-il.Yl  ^(bz)]  .  (D.14) 
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With  the  property 


2 

{f-[Kv2Uz)  -  K^UzjK^Uz)]}’  -  zKv2(Cz)  , 


the  normalization  integral  is  carried  out  to  give 


N  -  —  .  *  -  [K  ,(bz.)K  .(bz.)]'*  . 

b  zQ  2sinhp  v-i  0  v+1  0 


(D. 15) 


(D. 16) 


It  is  noted  that  kv(bZg)  »  0  due  to  (D.9)  and  (D.14).  With  the  use  of  (D.14) 
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and  the  recursion  formula 

Iv-l(*)  -  Iv+i^z)  *  (2v/z)Iv(z)  ,  (D. 17) 


(D.16)  can  be  written  as 


(D.  18 


Nb  ■  ^WV  *  Iv-l(bzO)l 


-1 


which  is  more  suitable  for  numerical  calculations  based  on  (D.ll).  In 
particular,  in  the  limit  b  •+  0,  from  (D.ll),  (D.12)  and  (D.18),  we  have 


Nfa  «  b|l  -  (bzQ/2)2/ v(  1  -v)  |  ' 1  <*  b  . 


(D. 19; 


Appendix  E.  Evaluation  of  (4.18) 

Here  we  give  some  technical  remarks  on  the  numerical  evaluation  of  (4.18). 
Using  E  =  kZg  and  E  =  bz^,  one  can  write  (4.18)  as 


M 


t 


(E)0(6,E)  +  c.c.  ]  , 


(E.l) 


where 


0(6, E) 


[ly(Ez)  -  I  ^(Ez))Jvj(Ez) 


(E.2) 


The  first  thing  to  do  is  to  replace  «  by  a  finite  z^  in  (E.2).  For  a  given  E, 
z^  is  roughly  a  classical  turning  point  of  the  potential  (4.11), 

zm  s  /y/E  .  (E. 3) 

m 


Next  we  evaluate  the  Bessel  functions  I  *1*  and  J  for  the  variable  ranges  E 

v  -v  v  ° 

<  Ex  <  /y  and  E  <  Ez  <  Ez  .  For  a  realistic  situation  /y  <  10,  the  Bessel 

—  —  n  - 

function  1^  can  -a  evaluated  accurately  by  the  power  series  expansion  (D.ll). 

On  the  other  hand,  even  for  an  interesting  situation,  E  -  1,  Ez  a  E/y/E  can  be 

ID 

large  for  small  6.  We  have  found  that  the  power  series  expansion  (D.2)  is  good 

18 

for  z  <  30.  For  z  >  30,  one  can  use  Hankel's  asymptotic  expansion, 

Jvu>  c°s(z  -  — —  ~ 1  If)  -  Q(v,z)sin(z  -  n))  ,  (E.4) 

2 

where  with  £  =  4v  , 


Though  it  is  a  bit  time-consuming,  the  integral  representation 


18 


IvU)  = 


js/zr, 

/■rr r ( v  f  j) 


■f 

J0 


d8  cos(zcos8) 


•  2vo 
sin  0 


(E.7) 


is  good  for  an  arbitrary  argument  z.  A  caution  for  the  integration  in  (E.7)  is 
that  at  small  0  <<  1 ,  the  integrand  oscillates  with  a  period  -  0. 
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Figure  Captions 

1.  (a)  Three-dimensional  geometry  for  the  scattering  eigenstate  characterized 

by  the  parallel  and  perpendicular  wave  numbers,  k  and  k  . 

X  z 

(b)  One-dimensional  simplification  when  the  parallel  and  perpendicular 
motions  are  approximately  separable. 

(c)  Two-dimensional  T-shape  model  for  low- temperature  adsorption. 

2.  Static  square-well  potential. 

3.  Schematic  figure  of  the  reduced  potential  U(x)  and  its  modifications  (dashed 
lines)  in  the  region  R*  <  x  <  R. 


4.  k£  as  a  function  of  the  exponent  n  for  y  *  20  and  g  •  0.6. 

5.  Self-trapping  model  for  He  adsorption  on  W. 

6.  One-particle  energy  spectra  for  the  surface  W-atom  and  He-atom. 

7.  (a)  | S |  and  ReS  vs.  k.  It  is  seen  that  S  +  1  as  k  +  0. 

(b)  a(k)  vs.  k.  Below  kc  -  0.1  A  a  has  an  approximate  linear  dependence 
on  k. 

t  *2 

8.  The  solid  curve  is  the  long-range  potential  (4.11)  in  the  unit  of  A  with 
Zq  *  0.5  A  and  p  *  5.0.  The  horizontal  solid  lines  denote  its  bound  states 
calculated  from  (D.9).  Bound  states  with  very  small  energies  are  described 
by  (D. 12) . 

9.  a(k)  vs.  k  for  D  *  1,  2  and  3.  The  parameters  are  Zq  *  0.5  A  and  p  «  5.0. 
For  the  D  *  1  case,  a  spurious  logarithmic  singularity  is  seen  for  small  k, 
as  described  in  the  text.  The  other  cases  D  *  2  and  3  demonstrate  finite 
quantum  sticking  coefficients  at  k  -*■  0. 
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